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The origin of large scale magnetic fields in astrophysical rotators, and the conversion of gravitational energy into radiation 
near stars and compact objects via accretion have been subjects of active research for a half century. Magnetohydrodynamic 
turbulence makes both problems highly nonlinear, so both subjects have benefitted from numerical simulations. However, 
understanding the key principles and practical modeling of observations warrants testable semi-analytic mean field theo- 
ries that distill the essential physics. Mean field dynamo (MFD) theory and alpha-viscosity accretion disc theory exemplify 
this pursuit. That the latter is a mean field theory is not always made explicit but the combination of turbulence and global 
symmetry imply such. The more commonly explicit presentation of assumptions in 20th century textbook MFDT has 
exposed it to arguably more widespread criticism than incurred by 20th century alpha-accretion theory despite comple- 
mentary weaknesses. In the 21st century however, MFDT has experienced a breakthrough with a dynamical saturation 
theory that consistently agrees with simulations. Such has not yet occurred in accretion disc theory, though progress is 
emerging. Ironically however, for accretion engines, MFDT and accretion theory are presently two artificially uncoupled 
pieces of what should be a single coupled theory. Large scale fields and accretion flows are dynamically intertwined be- 
cause large scale fields likely play a key role in angular momentum transport. I discuss and synthesize aspects of recent 
progress in MFDT and accretion disc theory to suggest why the two likely conspire in a unified theory. 
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1 Introduction 

Large scale dynamo (LSD) theory in astrophysics is aimed 
at quantitatively understanding the in situ physics of mag- 
netic field growth, saturation, and sustenance on time or spa- 
tial scales large compared to the turbulent scales of the host. 
The presence of magnetohydrodynamic turbulence makes 
the problem highly nonlinear, exacerbating the importance 
of numerical simulations. However, large scale spatial sym- 
metry and slow evolution of large scale fields compared to 
turbulent fluctuations time scales motivates a mean field ap- 
proach in which statistical, spatial, or temporal averages 
are taken and evolution of the mean field studied (Mof- 
fatt 1978). An important goal of 21st century semi-analytic 
mean field dynamo theory (MFDT) is to capture the essen- 
tial nonlinear physics in a minimalist way that is simple 
enough for phenomenological models, but also accounts for 
the nonlinear saturation found in simulations (for technical 
reviews see: Brandenburg & Subramanian 2005; Blackman 
2007). The past decade has revealed significant progress in 
this endeavor as the growth and saturation of LSDs seen in 
simulations seems to be reasonably matched by mean field 
dynamo (MFD) theories that follow the time dependent dy- 
namical evolution of magnetic helicity (e.g. Blackman & 
Field 2002). 
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Stars show that in situ LSDs must operate in nature. 
The sun exhibits sign reversals of the large scale field over 
each half solar cycle (e.g. Wang & Sheely 1993; Schrijver & 
Zwaan 2000). If the field were simply that frozen into to the 
ISM during initial formation, there would be no reversals. 
LSDs are also likely operating in galaxies, despite whatever 
initial fields may have been seeded cosmologically because 
supernova driven turbulence requires sustenance of the large 
field against turbulent diffusion. Turbulent diffusion in 3- 
D is not likely suppressed as it may be in in 2-D simula- 
tions for weak magnetic fields (Cattaneo 1994) and there is 
little evidence that turbulent diffusion in interstellar media 
is quenched. Detailed modeling of the geometry of large 
scale fields can be matched by mean field theories (Beck 
et al. 1996,Vallee 2004, Shukurov 2002), though incorpo- 
ration of the nonlinear LSD principles of the past decade 
for idealized dynamos into more geometrically realistic dy- 
namos of astrophysical rotators with realistic boundary con- 
ditions setting is an ongoing avenue of research (e.g. Sur et 
al. 2007,2009). 

Given the role of LSDs in the turbulent rotators of stars 
and galaxies, their presence in turbulent accretion discs is 
also likely. The large scale fields required by prevailing mod- 
els of astrophysical jets may be produced by an in situ dy- 
namo, although understanding the relative importance of in 
situ production vs. flux accretion in that context has been a 
long-standing topic of research (e.g. Lovelace et al. 2009). 
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As in stellar and galactic environments of LSDs, the likely 
presence of magnetohydrodynamic turbulence in accretion 
discs makes the problem fundamentally nonlinear and a fo- 
cus of numerical simulations. Computational resources have 
demanded local shearing boxes for thin disc studies but, 
as discussed in section 4, even localized stratified shear- 
ing box simulations show evidence for LSDs and the in- 
fluence of large scale fields (Brandenburg et al. 1995; Lesur 
& Ogilive 2008; Davis et al. 2009). Non-thermal coronal 
luminosity from accretion engines also hints at the role of 
large scale field because larger scale structures more eas- 
ily survive the buoyant rise before dissipating in coronae 
(Blackman & Pessah 2009). 

While LSD theory has always been presented with ex- 
plicit mean field theories that expose the approximations 
made, axisymmetric accretion disc models such as the Shakura 
Sunyaev (Shakura & Syunyaev(1973), hereafter SS73) vis- 
cous a ss model are not always recognized as mean field 
theories by those using them to model spectra. In this sense, 
21st century LSD theory is more progressed than axisym- 
metric mean field accretion disc theory: The a ss disc pre- 
scription does not yet offer predictive power for numerical 
simulations, whereas 21st century LSD theory predicts the 
saturation of LSDs seen in a range of simulations. Ironically, 
because large scale fields likely play a fundamental role in 
accretion discs, both MFD and accretion theory are really 
two faces of a single coupled theory. 

Herein, I first explore some conceptual progress and un- 
derlying principles of 2 1 st century LSD/MFD theory in broad 
brush strokes. I then discuss parallels between large scale 
dynamo theory and accretion disc theory as separate en- 
deavors. Finally, I bring the two together and identify evolv- 
ing connections and the importance of integrating the two 
theories into one. 

2 Overview of 21st Century LSD Theory 

For ~ 50 years, 20th century textbook MFD theory (e.g. 
Moffatt (1978)) has lacked a saturation theory to predict 
how strong the large scale fields get before quenching via 
the back-reaction of the field on the driving flow. But sub- 
stantial progress toward a nonlinear mean field theory has 
emerged in the 21st century via a symbiosis between analyt- 
ical and numerical work. Coupling the dynamical evolution 
of magnetic helicity into the dynamo equations turns out to 
provide fundamental insight to the saturation seen in simu- 
lations. Whether this approach is the only way to understand 
the saturation is uncertain but it has been successful. 

The connection between magnetic helicity and large scale 
dynamos is implicit in the classic EDQNM spectral model 
of helical MHD turbulence in Pouquet et al . ( 1 976) . In Klee- 
orin & Ruzmaikin (1982) an equation that, in hindsight, 
couples the time evolution of mean magnetic helicity to the 
dynamo is present though not actually solved for the time 
evolution. The success of the coupling becomes dramati- 
cally evident only when the time evolution is studied. That 
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exploration has led to the paradigm shift in MFD theory. For 
a sampling of such papers see Blackman & Field (2002); 
Brandenburg & Subramanian (2005); Blackman (2007); Kapyla 
et al. 2008). 

2.1 Globally Helical 

Most work in LSD/MFD theory has focused on systems 
that are initially globally reflection asymmetric (GRA). In 
such systems, pseudoscalars, such as the hemispherically 
averaged product of angular velocity and density gradient 
imposed by the initial setup facilitate e.g. a kinetic helic- 
ity pseudoscalar, common to standard textbooks (Moffatt 
(1978),Parker 1979 ) "ad yn effect" of mean field dynamos. 
However, 21st century MFD theory has revealed that it is 
the magnetic helicity which should be considered the unify- 
ing quantity for all MFDs. The initial GRA conditions and 
the potential pseudoscalars that arise are important primar- 
ily because they provide one way to sustain an electromotive 
force aligned with the mean magnetic field. This alignment 
is a source of magnetic helicity or magnetic helicity flux. 

There are two classes of GRA LSDs. The first is flow 
driven helical dynamos (FDHD) which apply inside of as- 
trophysical rotators. Here the initial energy is dominated by 
velocity flows and the field responds. These are potentially 
linked to coronae where a second type of LSD, the mag- 
netically driven helical dynamo (MDHD) can operate. A 
MDHD can also be described as dynamical magnetic relax- 
ation (Blackman & Field 2004) toward the "Taylor" state 
(Taylor 1986). Note that the Taylor state is the lowest en- 
ergy state to which a magnetically dominated configuration 
evolves. For a closed system, the state is characterized by 
the magnetic helicity migrating to the largest scale avail- 
able subject to boundary conditions. In an accretion disc, 
the MDHD would characterize coronal relaxation of mag- 
netic structures (fed by helical fields from below) into larger 
(even jet mediating) scales in a magnetically dominated en- 
vironment. In this respect, the MDHD in disc-corona or star- 
corona interfaces are analogous to laboratory plasma dy- 
namos (Bellan 2000, Ji & Prager 2002) that occur in reverse 
field pinches (RFPs) and Spheromaks with the corona being 
the internal volume of these configurations and the disc or 
stellar surface being the boundary. 

LSDs always involve some helical growth of the large 
scale field which is coupled to a helical growth of smaller 
scale fields of opposite sign and/or a compensating helic- 
ity flux when boundary conditions allow. When small scale 
magnetic or current helicity evolution is coupled to the large 
scale field growth, the simplest 21st century dynamo reveals 
that the mean field 'dynamo ctdyn becomes the difference 
between kinetic helicity and current helicity: For an a 2 dyn 
type FDHD simulated in a closed box (Brandenburg 2001), 
the current helicity builds up as the large scale field grows 
and quenches the FDHD, in accordance with MFD predic- 
tions (Blackman & Field 2002). 

The effect of boundary or flux terms can vary depending 
on the sign and relative flux of small and large scale mag- 
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netic helicities. For a closed system the buildup of small 
scale magnetic (and current) helicity is necessarily accom- 
panied by small scale magnetic helicity buildup with the op- 
posite sign, which quenches the LSD to at most a rate deter- 
mined by the dreadfully slow resistive dissipation of small- 
scale helicity. If this "catastrophic" quenching occurs before 
enough large scale field is grown in the fast growth regime, 
a preferential flux of helicity of the opposite sign through 
a realistic astrophysical boundary is desirable to alleviate 
this quenching and sustain further fast fast growth. For GRA 
systems, in each hemisphere, the large scale helicity builds 
with one sign and the small scale with the opposite sign so 
the flux of the small scale helicity is desired. (Blackman & 
Field 2000; Vishniac Cho 2001; Sur et al. 2007; Kiipyla et 
al 2008). It is possible that even if the field grows to accept- 
able magnitudes without boundary fluxes, the cycle period 
becomes resistively limited. In this case the helicity fluxes 
can unclog the cycle. 

The actual mechanism by which small scale helicity might 
be preferentially ejected (rather than larger scale magnetic 
helicity) depends on the astrophysics of particular system. 
The physics of helicity ejection should emerge as part of 
the natural buoyant loss of magnetic flux from a rotator, and 
accordingly in simulations. Numerical simulations support 
these basic principles: when shear is present, and when sur- 
faces of constant shear align toward open boundaries a spe- 
cific non-vanishing magnetic helicity flux is allowed which 
alleviates the quenching (Brandenburg & Sandin 2004; Kapyla 
et al. 2008). 

For an MDHD, the system is first dominated instead by 
the current helicity and a growing kinetic helicity can act 
as the back-reaction (Blackman & Field 2004). Both sim- 
ple FDHDs and MDHDs are accessible within the same 
formalism, all unified by tracking magnetic helicity evo- 
lution, and aided by thinking of the field as ribbons rather 
than lines (Blackman & Brandenburg 2003). More work on 
how the fields evolve from within the rotator to produce the 
global scale fields in coronae is needed. Interestingly, al- 
though most observations probe coronare, most work on as- 
trophysical dynamo theory has focused on the FDHD and 
not the MDHD. 

2.2 Not Globally Helical 

LSD action has also been observed in non GRA simula- 
tions (e.g. Yousef et al 08; Lesur & Ogilvie 08; Branden- 
burg et al. 2008). The minimum global ingredients for this 
class of LSD seem to be shear plus turbulence (Yousef et 
al. 2008). The non GRA LSDs grow large scale fields on 
scales larger than the outer scale of turbulence but smaller 
than the global scale. Although there is no GRA, in regions 
where the large scale field is coherent regions, there is a field 
aligned electromotive force (EMF), and thus an intermedi- 
ate scale source of magnetic helicity that may switch signs 
between coherence regions and globally averages to zero. 
It may be that the non-GRA LSD action always involves a 
local helicity flux between coherence regions. It is therefore 



important to recognize that the absence of GRA need NOT 
imply the irrelevance of magnetic helicity and its flux. 

In the absence of invoking GRA, a number of approaches 
that appeal to the anisotropy of the turbulence due to shear 
have been studied, and can provide finite contributions to 
the EMF on intermediate scales (e.g. Vishniac & Branden- 
burg 1997 [stochastic a^ yn effect]; Blackman 1998a [turbu- 
lent small scale dynamo + large scale shear] shear); Klee- 
orin & Rogachevskii 2003 [shear current effect, though see 
Brandenburg et al 2008 and particularly Sridhar and Subra- 
manian 2009ab which challenge its existence] Schekochi- 
hin et al 2008 [turbulence + shear]). The role of some kind 
of helicity flux (e.g. Vishniac & Cho 2001) between quasi- 
local sectors rather than globally is likely important, but the 
dominant mechanism and the consequences of the type of 
mean fields produced by non-GRA LSDs remain topics of 
desired research. 

3 Why Magnetic Helicity Is Relevant 

Magnetic helicity has emerged as a key quantity in under- 
standing LSD saturation but why should this be? Several 
concepts conspire to provide insight into why. 

3.1 Standard textbook mean field dynamos 
necessarily generate large scale magnetic helicity 

The mean field B satisfies the induction equation (Moffatt 
1978; Krause & Radler 1980) 

dB - 

¥ = - cVxE ' (1) 

where 

cE = -V x B - (v x b) + v M V x B, (2) 

2 

where the magnetic diffusivity vm = ^ and rj is the resis- 
tivity. The turbulent electromotive force £ = (v x b) can 
in general be expanded in terms of powers of spatial and 
temporal derivates of the mean magnetic field and mean ve- 
locity (or mean shear). The lowest order contributions to 
£ from the mean magnetic field terms lead to exponential 
growth and diffusion in 20th century textbooks (and are the 
only contributions in the absence of a mean velocity) are 
given by 

£i = onjBj — PijkdjBk- (3) 

Using Maxwell's equations, the definitions VxA = B, 
and J = ^ V x B, along with vector identities, the equation 
for the evolution of magnetic helicity density H M = ( A-B) 
is (see Bellan 2000; Blackman 2007) 

dt(A-B) = -2i/ A/ (J-B)-V-(2$B + Ax9 t A),(4) 

where $ satisfies E = V<I> — -%r- Following the analogous 
procedure for the mean and fluctuating quantities, and using 
(|2]) and (J3]) gives for the large and small scale quantities 
Hf 4 = (A ■ B) and = (a ■ b) respectively 

d t (A-B) = 25-B-2z/ M J-B-V-(2$B+AxdtA)(5) 
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and 



d t a b = -2£-B-2z/ M j ■ b- V • (2^b + a x d t a).(6) 

To highlight that standard large scale dynamos neces- 
sarily involve generation of mean magnetic helicity, note 
that the minimalist versions of such dynamos invoke Oij — 
and fa), — v M,T^ijk where adyn is a pseudoscalar 
and vm.t acts as a mean field scalar diffusion coefficient. 
(Note that in anisotropic systems, /Jy^ may have negative 
components that might act to amplify rather than diffuse-an 
ingredient in some of the models of non GRA LSDs dis- 
cussed in section 2.2.) The growth of the large scale field 
then involves a finite S ■ B, which provides an equal and 
opposite source term to the large scale and small scale mag- 
netic helicity as seen from (|5]l and |6]). In solving stan- 
dard 20th century textbook theory focuses on an equation 
for the large scale field from which one can write an equa- 
tion like |5]). But 20th century dynamo theory does not cou- 
ple in the evolution of the small scale magnetic helicity con- 
tained in {5). 

In the simplest of 21st century dynamos (e.g. Blackman 

& Field 2002), a dyn r«v • Vxv) - (b • Vxb» = 

OLkin + o-mag, where r is of order a correlation time of dom- 
inant fluctuations. In 20th century theory only a^in con- 
tributes to adyn- The contribution a mag is proportional to 
the small scale current helicity, which in turn is proportional 
to i?2 1 m tne Coulomb gauge for a closed system. There- 
fore the time evolution of ad yn is directly coupled to the 
time evolution of magnetic helicity. If Hi grows initially 
from a^in, then H^ 1 grows of opposite sign which in turn 
quenches adyn- As alluded to in our discussion of MDHDs 
above, it is also possible to drive the dynamo instead with 
ctmag, such that akin is the quencher. I come back to this 
point in section 5. 

A subtle aspect of the helicity density equations |5]) and 
|6]l is the issue of gauge invariance. In the absence of bound- 
ary flux terms, all terms are manifestly gauge invariant. In 
the present of flux terms, the gauge invariance is broken 
(e.g. Berger & Field 1986). Subramanian & Brandenburg 
(2006) carefully construct a generalized local helicity den- 
sity that reduces to the above equations in the absence of 
flux terms and has the same form in the presence of flux 
terms but with a redefinition of Hf 1 , H^ 1 and the flux terms 
that avoids the use of the vector potential. The lowest order 
terms in their definition are similar to what is obtained from 
the Coulomb gauge but the corrections make each individ- 
ual terms in the equation gauge invariant with a local phys- 
ical meaning. That important construction puts the physical 
discussion of helicity density evolution equations here and 
elsewhere on firm ground, particularly if one seeks to mea- 
sure a local flux density as a physical quantity. However, 
if one uses the helicity evolution equations as intermediary 
equations in a particular gauge and can convert back to the 
magnetic field, the resulting magnetic field will always be 
gauge independent. 



3.2 Evolution of Total Magnetic Helicity vs. Total 
Magnetic Energy 

The Taylor state (Taylor 1986), mentioned in section 2.1, 
is typically derived by starting with a bounded, static mag- 
netically dominated configuration and minimizing the mag- 
netic energy subject to the constraint that magnetic helicity 
is conserved. The result is a configuration in which the mag- 
netic helicity evolves to the largest scale available subject to 
boundary conditions. This is NOT a dynamical calculation, 
but a calculation of an end state. In fact, the evolution to 
this end state is a dynamo-albeit a MDHD. Recall that the 
MDHD is an LSD for which magnetic helicity is injected 
at small scales and (non-locally) inverse cascades to large 
scales (Bellan 2000; Blackman & Field 2004). 

The concept of the Taylor state hints at a fundamental 
feature of magnetic helicity and its connection to large scale 
dynamos. Even a FDHD in a closed volume involves an in- 
verse transfer of magnetic helicity to large scales (Pouquet 
et al. 1976; Brandenburg 2001) For a system with initially 
zero net magnetic helicity the FDHD segregates one sign 
to large scales (e.g. the box scale) and the other sign to the 
largest scale available to it (the forcing scale). The system 
then achieves a "doubly relaxed state" (Blackman 2003). 

Key in deriving Taylor states is the assumption that mag- 
netic helicity evolves slowly compared to magnetic energy. 
This also emerges as justifiable a posteriori for numerical 
simulations of the FDHD (e.g. Brandenburg 2001). It is there- 
fore instructive to evaluate how widely applicable this as- 
sumption can be expected to hold. More specifically, one 
can ask for what range of magnetic energy and magnetic 
helicity spectra does the the magnetic helicity decay more 
slowly than magnetic energy? 

Following Blackman (2004) and working in the Coulomb 
gauge, we write the total magnetic energy density 



(B ) =M = 



M k dk 



(7) 



where k and k Vu are the minimum and maximum (resis- 
tive) wave numbers, and where the magnetic energy density 
spectrum 



\B\ 2 k 2 dn k 



A| 2 /c 4 dfL oc fc- 



Here fij. is the solid angle in wave-number space, q is as- 
sumed constant, and the tilde indicate Fourier transforms. 
The magnetic helicity density spectrum is then 



H k = \j A(fc)B*(fc) + A*(fc)B(fc) 



k 2 dQ,k 



M k f(k)/k, 



(9) 



where f(k) oc k~ s is the fraction of magnetic energy that 
is helical at each wave number and s is taken as a constant. 
We then also have correspondingly 

(A-B)= / /(^Mkk^dk (10) 
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and the current helicity 



(J B) 



f(k)kM k dk. 



(11) 



Now, using ([5]) for a closed system we have 

8tH = dt(A'B) = -2v M (J-B) (12) 

and the magnetic induction equation gives 

d t M = d t (B 2 ) = -2 l / M ((VB) 2 )-2(v(JxB)>.(13) 

Using ([10]), (fTTj, ([12) and ([13} we then obtain 



th 



-H 



f k l" f{k)M k k~ l dk 
2 "m jt M f(k)kM k dk 



and 



—M 



tm = 



kT M k dk 



(14) 



(15) 



(d t M) res 2v M $t L M k 2 M k dk ' 
where the subscript "res" indicates the contribution from 



the penultimate term in ( 13 i only. The range of s and q 
for which R = ^ > 1 corresponds to regime in which 
the magnetic helicity decays more slowly than the magnetic 
energy. Blackman (2004) showed that R > 1 for the com- 
bination of s > and 3 > q > 0. and that R < 1 for 
small < q < 1 and s < 0. Standard derivations of the 
Taylor state are correct only when R > 1 . Indeed R > 1 for 
a wide range of reasonable spectra and magnetic helicity is 
therefore often conserved better than magnetic energy. 

Two key implications of this last statement are that: (1) 
the evolution of magnetic helicity in a system, its separation 
of scales, its flow out of a boundary, and its rate of dissi- 
pation should be tracked in dynamo models and is missing 
from 20th century textbook dynamos which do not conserve 
helicity (see Fig 1 . of Blackman & Brandenburg 2003); and 
(2) Magnetic helicity conservation should really be consid- 
ered a more fundamental constraint on closed box numerical 
MHD simulations than magnetic energy conservation. This 
also applies to shearing box simulations (section 4). 



4 Accretion Theory: Role of Mean Field 
Theory and Large Scale Fields 

4.1 Standard accretion theory is a mean field theory 

Gaseous accretion discs around astrophysical rotators have 
long been studied as a source of luminosity for compact 
engines (see Treves et al. 1989; Frank et al. 2002). In as- 
tronomy, a practical semi-analytic formalism to be incor- 
porated into emission and spectral models for comparison 
with observations is a desired goal. Because of this goal, 
the presentation of practical accretion models in standard 
texts (e.g. Frank et al. 2002) based on SS73 employs the 
Navier-Stokes equation and replaces the microphysical vis- 
cosity with a turbulent viscosity without explaining the the- 
oretical complexity and assumptions underpinning this bold 
procedure. The replacement of the microphysical viscosity 



by a turbulent viscosity is itself a turbulent closure and war- 
rants a formal derivation analogous to that of the MFD the- 
ory coefficients which also require closures. Closures ap- 
proximate an otherwise infinite set of nonlinear turbulence 
equations by an approximation that facilitates a finite set. 

That the SS73 model is both axisymmetric AND in- 
volves a turbulent viscosity implies that the theory is a mean 
field theory. Local axisymmetry never applies on the scale 
of turbulent eddies. The sense in which accretion disc the- 
ory is a mean field theory is addressed in a subset of theoret- 
ically oriented literature (e.g. Balbus et al. 1994; Blackman 
1997; Balbus & Hawley 1998; Ogilvie 2003 Pessah et al. 
2006; Hubbard & Blackman 2009). In contrast, 20th cen- 
tury textbook LSD theory is typically more explicitly pre- 
sented as a mean field theory (e.g. Moffatt 1978, Parker 
1979). This has arguably led to the 20th MFD theory being 
subject to broader critical scrutiny than SS73 despite sim- 
ilar weaknesses. MFD theory has traditionally focused on 
the magnetic induction equation as the "primary equation" 
of the theory and accretion disc models focus on the mass 
conservation and angular momentum transport equations as 
the "primary equations" of the theory with the coupling to 
the magnetic field often ignored or "swept" into "a ss ". 

A practical consequence of accretion theory being a mean 
field theory is a built in limitation in predictive precision as- 
sociated with the stochasticity (Blackman 1998b). For ob- 
servations that collect data on time scales shorter than the 
largest eddy turnover time divided by the square root of the 
number of eddy cells contained in a spatially unresolved re- 
gion contributing emission within a given frequency range, 
a mean field theory does not make an exact prediction. There 
is an intrinsic variability associated with predictions from 
any mean field theory. One has to be sure that the interpreta- 
tion and use the theory corresponds with its correct meaning 
when compared to a given set of observations. 

At present, the gap between formal identification of the 
minimal set of transport coefficients that best captures the 
angular momentum transport from turbulence in accretion 
theory seen in simulation and practical models for use in 
spectral modeling, is larger than the analogous gap in MFD 
theory. MFD theory now has a predictive semi-analytic the- 
ory that agrees with a set of nonlinear simulations which 
is not the case for mean field accretion theory like SS73. 
MFD theory has presently 2 tensor coefficients in common 
use (the and fiu k coefficients which are used even in 
semi-analytic dynamo models to match observed field struc- 
tures in stars or in galaxies. Dynamical equations for these 
quantities that incorporate magnetic helicity evolution and 
the "minimal r" closure (e.g. Blackman & Field 2002, see 
also Snellman et al. 2009 applied to shear flows) capture 
MFD saturation. In accretion theory, there are few models 
that use anything other than the single scalar diffusion coef- 
ficient of the SS73 formalism. The minimum properties that 
turbulence must have to transport angular momentum out- 
ward are more subtle than what purely isotropic turbulence 
provides, and it is likely that more than one transport coef- 
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ficient is needed (Pessah et al. 2006; Hubbard & Blackman 
2009). Progress toward capturing the saturated stress of the 
magneto-rotational instability (MRI) as a mean field model 
using a closure similar to the "minimal r" approximation 
and going beyond the a ss disc model is emerging (Ogilvie 
2003; Pessah et al. 2006). 

Although the magneto-rotational instability (MRI) has 
emerged as a leading candidate for local angular momentum 
transport in accretion discs (e.g. Balbus & Hawley 1998; 
Balbus 2003), there remains a disconnect between what shear- 
ing box simulations have told us about the MRI vs. how the 
instability might operate in nature. To date, simulations have 
told us that the MRI plausibly operates in accretion disc 
but not robust scalings of transport coefficients that spec- 
tral modelers can use. Patience (likely for for several more 
decades) is required, as the computational and conceptual 
demands are substantial. 

To achieve maximal resolution, most first generation MRI 
simulations (except Brandenburg et al.1995) did not use ex- 
plicit viscosity or magnetic diffusivity. Indeed recent work 
does show that for unstratified simulations there is signifi- 
cant magnetic Prandtl number dependence (Fromang et al. 
2007), though this dependence is reduced for stratified boxes 
(Davis et al. 2009; Shi et al. 2009). Similarly, the angu- 
lar momentum transport coefficient a ss depends strongly 
on the box size and the strength of the initially imposed 
weak mean field strength for unstratified boxes with im- 
posed mean fields (Pessah et al. 2007). Interestingly, a ss 
varies <~ 4 orders between simulations, but a S8 /3 p , where 
(3 P is the ratio of thermal to magnetic pressure, is nearly 
a constant (Blackman et al. 2007). A theory that explains 
how the MRI saturation depends upon all these choices is 
lacking, thereby limiting our quantitative implications for 
understanding angular momentum transport in a realistic 
system. For unstratified boxes, the dependence of the MRI 
saturation on parasitic modes and the need to have large 
enough boxes to accommodate them is possibly quite im- 
portant (Pessah & Goodman 2009). The mean field approach 
of Pessah et al. (2006) hints at what minimal set of coeffi- 
cients might be required in a more robust but still practical 
model of MRI mediated accretion discs. The approach bet- 
ter matches the stress dependence on rotation profile than 
the a ss of SS73 in numerical simulations of unstratified 
boxes with an initial mean field (Pessah et al. 2008). . 

4.2 Why the focus on local transport? 

Accretion disc theory highlights is a duality of influence of 
computational tools. 

The practicality of the SS73 model has focused atten- 
tion toward local models of angular momentum transport. 
Indeed, for much of the past 15 years, the promising MRI 
has been studied as a local instability or "source" of the local 
turbulence that can transport angular momentum outward. 
The computational focus has employed shearing boxes where 
a small radial slab of the disc is considered (Balbus & Haw- 



ley 1998, Balbus 2003). The focus on local boxes is also ex- 
acerbated by computational limitations. It has not yet been 
possible to simulate a global thin disc. 

Understanding the local physics is important but it is 
also noteworthy that some paradigms of accretion discs in- 
volved buoyant magnetic loops (Lynden-Bell 1969; Field 
& Rogers 1993; Tout & Pringle 1992; Johansen and Levin 
2008) or outflows (Blandford & Payne 1982; Konigl 1989, 
even Colgate et al. 2001) involving large global scale mag- 
netic fields that transport the angular momentum. Maybe 
these are the primary modes of transport. 

In short, although numerical simulations have become 
fundamental tools for progress, the limitations of tools can 
also limit the focus. Large scale and nonlocal features and 
outflows may dominate the local transport. There has been 
some renewed study along these lines (Kuncic & Bicknell 
2007; Dobbie et al. 2008). Hints that the role of large scale 
structures may be quite fundamental is discussed next. 

4.3 Role of Large Scale Fields: Clues Even from 
quasi-local MRI Simulations 

Global simulations of thin accretion discs are a widely de- 
sired community goal. Based on the sun and the observation 
of jets from accretion engines and the existence of energet- 
ically significant coronae, large scale magnetic structures 
and outflows must be playing some role in non-local angular 
momentum transport in non-self gravitating discs. However 
even the study of quasi-local simulations is providing clues 
that non-local processes are important. 

First, note that recently a systematic study of magnetic 
Prandtl number dependence of MRI saturation has begun 
(e.g. Fromang et al. 2007; Davis et al. 2009). This is an on- 
going enterprise but noteworthy is that stratified simulations 
seem to show less dependence on its value than unstratified 
solutions (Davis et al. 2009). The stratified simulations also 
show less dependence on box size than do unstratified simu- 
lations (Pessah et al. 2007; Bodo et al. 2008). (How the role 
of parasitic modes as a possible mechanism of MRI satu- 
ration (Pessah & Goodman 2009) depends on box size for 
stratified vs. unstratified cases is uncertain.) 

Since real systems are stratified, the weaker dependence 
of the saturated state on box properties is step forward, but 
with stratification comes magnetic buoyancy and coronae. 
Large scale structures have the longest turbulent diffusion 
times and are most likely to survive the buoyant rise (Black- 
man & Pessah 2009). Thus if stratification is important then 
understanding what radial and vertical scales are large enough 
to fully capture the dominant magnetic stress remains a ques- 
tion that more global simulations may be needed to answer. 
Presently, the radial width of the Cartesian box simulations 
are significantly smaller than the disc radius but magnetic 
coupling of different regions in the disc via magnetic loops 
can non-locally transport angular momentum (e.g. Lynden- 
Bell 1969; Blandford & Payne 1982; Tout & Pringle 1992; 
Field and Rogers 1993). 



© 20 10 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journaI.org 



Astron. Nachr. / AN (2010) 



795 



Another point is that shearing box simulations show the 
generation of large scale mean toroidal magnetic fields that 
last many orbit times and exhibit sign reversals (Branden- 
burg et al. 1995, Lesur & Ogilvie 2008, Davis et al. 2009). 
These are seen in stratified simulations with and without pe- 
riodic vertical boundaries, and in unstratified simulations 
with periodic vertical boundaries (Lesur & Ogilvie 2008). 
The patterns indicate a large scale dynamo operating con- 
temporaneously with the small scale dynamo. When strat- 
ification is present, there are global (or hemispheric) pseu- 
doscalars which depend on the product of rotation and den- 
sity stratification that can supply a mean magnetic field aligned 
EMF (e.g. Brandenburg & Donner 1997; Brandenburg 1998; 
Tan & Blackman 2004) . But even in the unstratified cases, 
there seems to be a mean field aligned EMF in sub-global 
but non-local regions of the disc. Any presence of a mean 
field aligned EMF is then a sub-global but non-local source 
of magnetic helicity or magnetic helicity flux of opposite 
sign between sectors. All of this links back to the LSD dis- 
cussion of section 3. 

Finally, note that shearing box simulations impose a steady - 
state by artificially enforcing the rotation profile with no 
actual accretion. This rotation is not subject to the back- 
reaction of the amplified field or flow dynamics. Hubbard 
and Blackman (2009), argue that this may be too restrictive: 
In a real disc, energy in differential rotation is sustained only 
by accretion itself. If a steady-state is to be maintained via 
turbulent transport alone then there must be 100% power 
throughput from differential rotation to the turbulent cas- 
cade. This is not guaranteed if large scale fields drain power. 
A realistic steady-state solution would also have to then in- 
corporate stresses from large-scale fields. 

4.4 Coupling of Disc, Corona, and Jet 

That observations reveal jets and non-thermal emission best 
explained as emanating from coronae in accretion engines 
most directly highlights the importance of large scale fields 
therein. There is debate as to whether or not the very global 
scale fields of jets are accreted or produced in situ (Lovelace 
et al. 2009), but as emphasized, the sun, itself a stratified ro- 
tator, produces its global scale fields in situ because the sign 
of these fields reverse. Coronal loops are large scale with 
respect to the interior eddies but small scale with respect 
to the global scale fields of coronal holes. The formation 
of the global scale fields of coronal holes along which the 
solar wind propagates occurs as some of the closed loops 
open up (citesz). Jets from accretion engines in young stel- 
lar objects, active galactic nuclei and gamma-ray bursts or 
magnetic towers (Pudritz (2004), Lynden-Bell (2003), Uz- 
densky & MacFadyen (2006)) may be the direct analogue 
to coronal holes. These large scale coronal fields can be 
produced by the MDHD relaxation of smaller scale loops 
emerging from a FDHD inside the rotator below (Blackman 
2007) If the helical field injected from below is produced 
by a FDHD inside the rotator, the full description of the ori- 
gin of the large scale fields involves a coupling between the 



FDHD and MDHD through the coronal boundary, which in 
turn is coupled to the MRI in the accretion disc. 

Starting with the assumption that coronal luminosity from 
a turbulent accretion disc results from buoyant magnetic 
structures that survive turbulent shredding for at least one 
vertical density scale height, Blackman & Pessah (2009) de- 
rived lower limits on: (z) the scale of such magnetic struc- 
tures and (/;') the fraction of magnetic energy that needs to 
be produced above this scale l c within the disc to account for 
observed values of coronal to bolometric luminosity. They 
considered buoyant structures in pressure equilibrium with 
the ambient medium but with an additional magnetic en- 
ergy contribution from scales above l c , and thus a lower- 
than-ambient density. The results imply that typical coro- 
nal (non-thermal) to bolometric luminosity ratios observed 
in AGN require the critical scale for buoyancy to robustly 
exceed the scale of the largest turbulent motions l t and that 
double digit percentages of magnetic energy produced in the 
disc reside in large scale structures (e.g. flux tubes) whose 
smallest scale exceeds l t . 

The above results are consistent with recent work high- 
lighting the importance of in situ large scale dynamos in 
feeding coronae (Blackman & Field 2000; Vishniac 2009) 
and further motivate larger domains in stratified MRI simu- 
lations and determination of the magnetic energy and mag- 
netic helicity spectra produced therein. 

5 Combining MFD and Accretion Theory 
into a Single Mean Field Theory 

Understanding the competition between diffusion of global 
scale magnetic fields in accretion discs, advection of flux, 
and amplification of flux is all part of the same question: 
how do mean magnetic fields evolve in accretion discs when 
the dynamically correct mean field transport coefficients are 
obtained for both the coupled LSD and accretion dynamics? 
Ultimately, mean field accretion disc theory should be cou- 
pled to a MFD theory in a real disc as they are artificially 
separated components of what should be a single mean field 
theory. A coupled theory will involve separate the time evo- 
lution equations for mean velocity, surface density, mean 
magnetic field, and an energy equation. 

Some aspects of the coupled evolution of mean veloc- 
ity field and mean magnetic field in dynamo theory have 
been studied, though not part of accretion theory (Blackman 
1997; Blackman et al. 2006; Courvoiser et al. 2009). Ingre- 
dients of a more formal dynamical mean field accretion disc 
theory are also emerging (Pessah et al. 2007; Hubbard & 
Blackman 2009) which need to be combined with lessons 
learned form MFD theory. Campbell & Caunt (1999) com- 
bined 20th century MFD theory with the SS73 accretion 
theory closure to include local and large scale stresses. The 
SS73 type viscosity (via MRI) has also been used as a clo- 
sure for the turbulence from which semi analytic MFDT is 
then derived to produce large scale fields in discs to power 
outflows (e.g. Tan & Blackman 2004). 20th century MFD 
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theory has also been applied to MRI stratified simulations 
to explain observed cycle periods (Brandenburg et al. 1995; 
Brandenburg & Donner 1997; also Lesur & Ogilvie 2008; 
Davis et al. 2009). 

Although none of the above approaches yet combine 
insight gained from understanding nonlinear MFD theory 
via magnetic helicity conservation as discussed in section 
3 with accretion transport coefficients beyond that of SS73, 
there are indirect clues that the role of the a mag term dis- 
cussed in section 3.1 might be the driver of large scale field 
growth in sheared stratified rotators rather than the a^m 
term. This is based on (i) the sign of the a.d yn needed to 
model the parity and cycle periods seen in MRI simulations 
with a LSD theory (Brandenburg & Donner 1997; Branden- 
burg 1998, Rudiger & Pipin 2000; Rekowski et al. 2000; 
Gressel 2009), (ii) the mechanism of of a magnetic instabil- 
ity driven dynamo within stratified rotators (Spruit 2002), 
and (iii) the fact that the magnetic fluctuations seem to dom- 
inate the velocity fluctuations in all MRI simulations. With 
solutions to LSD models that match the accretion disk sim- 
ulations, the extent to which large scale fields transport an- 
gular momentum can be studied. 

The coupling of the disc to the corona in a real system 
highlights the role of magnetic energy flux and magnetic 
helicity fluxes so the dynamics of angular momentum trans- 
port in real discs are unlikely to be captured with unstratified 
simulations. 

6 Conclusions 

For any turbulent MHD system in which there is both turbu- 
lence on short time scales but ordered large scale patterns, 
there has to be a mean field theory that dynamically captures 
essential physics and key principles in a useful framework. 
Note that MHD itself is already a mean field theory which 
averages out the noise of individual particle motions. The 
question for turbulent rotators involving accretion and large 
scale dynamos is not "is mean field theory is correct?" but 
rather "do we have the correct mean field theory?" Progress 
in large scale dynamo theory is farther along that accretion 
theory in this regard. 

21st century mean field dynamo theory has incurred a 
paradigm shift compared to that of the 20th century by in- 
corporating the dynamical evolution of magnetic helicity. 
This has led to mean field models that successfully pre- 
dict the non-linear saturation of the large scale field growth 
seen in numerical simulations. In that way, 21st century dy- 
namo theory has overcome an long standing challenge of 
20th century dynamo theory. More work is still needed in 
applying the principles learned to realistic systems. 

The analogous paradigm shift has not yet occurred in ac- 
cretion theory. The single a ss viscosity parameter of SS73 
does not capture the quantitative properties of angular mo- 
mentum transport in numerical simulations or provide testable 
dynamical predictions for either local turbulent stresses pro- 
duced by the MRI, or large scale magnetic stresses. Some 



principles behind a generalized semi-analytic accretion model 
must include are emerging, and more effective closures than 
the a ss prescription are showing some promise. 

Ultimately, the need to incorporate both the large scale 
magnetic field growth and its stresses along with those pro- 
vided by small scale field, implies that semi-analytic mean 
field dynamo theory and semi-analytic accretion disc theory 
are two components of what is a unified mean field theory 
of accretion. Developing this unified theory is a challenging 
and meaningful goal for the present century. 
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